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Abstract. The presence of vortices in the ionisation of hydrogen atoms by positrons 
at intermediate impact energies is investigated. The present findings show that a 
previously reported minima in the fully-differential cross section is the signature of a 
vortex in the continuum positron-electron-proton system. The behaviour of the real 
and imaginary parts of the complex-valued transition matrix is studied in order to 
determine and characterize the vortex in momentum space. The obtained information 
is translated to fully-differential ionisation cross sections, feasible of being measured 
with currently available techniques. 



1. Introduction 

In the last few years, the study and understanding of positron induced reactions with 
atoms and molecules have benefited from a new breed of experiments (Surko et al 2005). 
The availability of intense mono-energetic positron beams and recent improvements in 
the measurement techniques are now routinely providing new evidence on such topics as 
Ps formation or positron annihilation and scattering. Besides the measurement of total 
and differential cross sections (e.g. Falke et al 1997, Marler and Surko 2005, Laricchia 
et al 2008), the future implementation of positron reaction microscopes would make it 
possible to obtain kinematically complete pictures of the product particles (Williams 
et al 2010, Mueller et al 2012). For instance, it might provide a way of analyzing 
such complex quantities as the fully differential cross section (FDCS) for positron- 
impact ionisation processes. It is noteworthy to point out that even in the simplest 
case of a hydrogenic target, the complete determination of the FDCS would require 
the measurement of four out of nine variables in the three-body momentum space, this 
latter number having been reduced by momentum and energy conservation and by the 
rotational symmetry about the collision axis. 

Even before any positron reaction microscope becomes operational, experimental 
and theoretical studies of one and two-body spectra have already leaded to a clear 
understanding of some of the main features of the FDCS for positron-impact ionisation 
collisions. This is the case, for instance, of the electron capture to the continuum (ECC) 
effect (Kover and Laricchia 1998), whose characteristics in the momentum spectra of 
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both the electron (Fiol and Barrachina 2011) and the recoil-ion (Barrachina and Fiol 
2012) are now well understood. 

However, up to our knowledge, the emergence of vortices in positron-impact 
ionisation collisions has not been researched neither experimentally nor theoretically. 
After being overlooked for decades, the presence of this kind of quantum structure in 
atomic and molecular processes was recently uncovered and even shown to be quite 
ubiquitous, being observed not only in electron-atom ionisation collisions (Macek et al 
2010, Colgan and Pindzola 2011), but also theoretically analyzed in the ionisation by 
the impact of ions (Macek et al 2009, Ovchinnikov et al 2011) and by intense electric 
pulses (Ovchinnikov et al 2010). However, as far as we know, vortices have not been 
predicted or even searched for in positron-impact ionisation collisions. 

It is noteworthy to mention that investigation of vortices in the case of electron 
impact presents additional complexities due to the indistinguishably of the projectile 
and target electrons. In the general case, the cross section may be written as the 
incoherent sum of the singlet and triplet terms. In the planar symmetric geometry 
studied in previous works one of these terms vanishes, and any vortex found in the 
other would produce a zero in the momentum distribution, as recognized earlier by 
Macek et al (2010). This is not longer true for the coUinear case chosen in our article 
but, on the other hand, neither is the explicit choice of a particular geometry necessary 
with positrons since there is not exchange contribution. 

More than two decades ago a deep minimum in the ionisation of atomic hydrogen by 
high-energy positrons was theoretically uncovered by Brauner and Briggs (1991). It was 
shown to appear when the electron and the positron emerged in the same direction at 45 
degrees from the forward direction. In principle, it was attributed to the interference of 
two indistinguishable processes akin to the double binary-collision proposed by Thomas 
in 1927. In the mechanism proposed by Brauner and Briggs (1991), the positron and 
the electron leave a first binary collision with the same speed, only to end up moving 
in a coUinear trajectory by the reorientation of any of them by an elastic scattering 
with the target nucleus. However, it is fair to mention that the minimum occurs at an 
electron speed smaller than that of the positron, an effect that increases for decreasing 
impact energies. 

At that time this structure was not recognized as a zero of the cross section, even 
though the sudden drop of about three orders of magnitude seemed to be pointing 
towards a zero rather than just a deep minimum. On the light of similar effects which 
had been experimentally (Murray and Read 1993) and theoretically (Berakdar and 
Briggs 1994) observed in (e,2e) collisions, and which only quite recently were recognized 
as vortices (Macek et al 2010), it can now be asserted with reasonable certainty that 
the effect theoretically discovered two decades ago by Brauner and Briggs (1991) in 
positron-impact ionisation collisions is a manifestation of a vortex. 

Unfortunately, since this effect was observed at 10 and 100 keV impact energies, 
and shown to be absent at 250 eV, it only seemed to occur at very large energies, where 
the FDCS will be too small to be measurable in the coincidence conditions required by 
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the collinear geometry (Brauner and Briggs 1991). 

Thus, some few years ago, we addressed the question of whether there was another 
unforeseen minimum in the FDCS that would persist at moderate impact energies and 
therefore would be amenable to experimental observation (Delia Picca et al 2004). 
Actually, there is a deep minimum located at precisely the energy and emission angle 
where an electron could emerge by lying in the saddle-point of the positron and the 
residual target-ion potentials. Unmistakable signatures of this structure were also 
pinpointed in the positron dispersion and the recoiling of the target ion (Delia Picca et 
al 2005). 

In the light of these evidences, the purpose of this article is to systematically 
investigate the precise location, structure and origin of this deep minimum and to 
demonstrate that it represents a fingerprint of a vortex in the transition matrix element 
(Navarrete et al 2011). 

2. Electron - Positron - Ion Three-Body Kinematics 

We investigate the emission of a single electron from an atom of ionisation energy e by 
the impact of a positron of initial velocity v Due to momentum conservation, the motion 
of the centcr-of-mass can be ignored and only two out of three positions or momenta are 
needed to completely describe the kinematics of the final three-body system. To this 
end, we will employ ion-centered coordinates, as defined by the relative positions r± 
and momenta k± of the electron (-) and the positron (-I-) with respect to the residual 
target of mass M. We will also use a set of Jacobi coordinates (Meyer 1999) given by 
the relative position r and momentum k of the electron with respect to the positron, 
and of their center-of-mass with respect to the residual target ion, namely R and K. 
The relation between the ion-centered and the Jacobi coordinates and momenta is given 
by r± = =F^/2 + R and k± ~ ^k + K /2, respectively. Here we are assuming that 
M ^ 1 (atomic units are used throughout this paper), so that r± and k± approximately 
represent the positions and momenta of the electron and the positron in a Laboratory 
reference frame. The fully differential cross section reads 



where Tij = (^j|Vi|^j) is the transition matrix clement. Here, the initial state 
accounts for the bound state of the target and the free motion of the positron with 
respect to the target. Vi is the interaction of the projectile with both particles in the 
target, namely Vi = Z/r_|_ — 1/r, where Z is the effective charge of the receding ion. Let 
us recall that the study of the three-body system with Coulomb interactions is one of the 
basic unsolved problems of Quantum Mechanics. Thus some kind of approximation is 
required for the description of the final continuum state ^J. Here wc choose to describe 
this final state by means of a correlated C3 wavefunction (Garibotti and Miraglia 1980, 




(1) 
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Brauner and Briggs 1986) 



(fc,r) 
^0 {k,r) 



(2) 



where 



= (27r)-3/2iV(;,/p) (-z^/p, 1; 



i(pq + p-g))e^ 



with normahzation Nlu) = r(l + iu) exp(7rz//2), is the continuum state for a two-body 
system of relative position q, momentum p and unity reduced mass, interacting via a 
Coulomb potential V{q) = — /i/q. 

Since the use of a C3 wavefunction is based on a perturbative approach it provides a 
better description of the ionisation process at large incident positron velocities. However, 
as explained in the introduction, the cross section magnitudes decrease strongly as the 
incident velocities increase, making very difficult any measurement of the FDCS. As a 
compromise, in this situation we will perform calculations at two collision energies: 100 
and 200 eV. While it is expected that the present theory would provide more accurate 
results at the larger energy, it is more likely that experimental confirmation will be 
reported at the lower energy of 100 eV. 

This approximation was already employed by us to successfully describe the electron 
capture to the continuum (ECC) effect in this same single ionisation collision by positron 
impact (Fiol et al 2001, Delia Picca et al 2006). Similar C3 approximations provided the 
first theoretical evidence of a deep minimum in positron-impact ionisation collisions at 
high energies (Brauner and Briggs 1991) and, with modified Sommerfeld u parameters, 
in (e,2e) processes (Berakdar and Briggs 1994). It also provided our first uncovering of a 
deep minimum of the FDCS at a saddle-point of the electron-positron-ion system in the 
ionisation of atoms by the impact of positrons of moderate energies (Delia Picca et al 
2004, 2005, 2006). Also, a variation of the C3 wavefunction, based on the works of Alt 
and Mukhamedzhanov (1993), modified to verify the correct asymptotic conditions has 
been employed to address (e,2e) processes (Berakdar and Briggs, 1994, 1994a), and very 
recently employed by Macek et al (2010) to investigate the appearance of vortices in 
electron-atom collisions in a symmetric geometry. It is noteworthy that this modification 
reduces to the C3 wavefunction (eq. [2]) for the collinear geometry employed in this work. 

3. Collinear Geometry 

Note that by energy conservation, the number of relevant scalar variables of the FDCS 
in Eq. ([T]) is reduced from six to five. Furthermore, in the absence of polarized targets or 
projectiles, the collision is axisymmetrical about the initial velocity v; and the dimension 
of the problem is further reduced to four. But even the mere representation of such a 
multidimensional object would be a formidable task, and this is so at two very basic 
levels. On one hand, there is the problem of choosing a set of four variables that would 
better represent the problem under study. Naturally, there is a large degree of freedom 
in the selection of these quantities. Some of these choices are related to different types 
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of experiments. This is the case of the so-called Rochester geometry (Baier et al 1969) 
where the cross section is written in terms of the polar angles of two outgoing particles, 
their relative azimuthal angle and the energy of one of them. This choice is complete in 
the sense that any other set of independent variables can be related to this one. 

The second difficulty with the FDCS is that it can not be represented in a 3-D 
space. The search for a practical solution to this conundrum is as old as Science and 
Art themselves (Malcolm 2004, Pauwels 2006). One way of making this representation 
manageable, but not the only option, is to reduce the dimension of the cross section. 
And there are two different ways of accomplishing this task. One option is to integrate 
the cross section over two or more of its variables. This is the case, for instance, in 
any single particle spectroscopy technique, where only the energy and/or the emission 
angle of one of the particles in the final state is investigated. The second option, and 
the one that we will adopt in this article, is to fix one or more variables to a particular 
value or condition. A very familiar and convenient choice when analyzing processes 
with two similar particles in the final state, as in (e, 2e) and (p, 2p) collisions (see, for 
instance, Lim and McCarthy 1966), is the so-called "symmetric geometry" introduced 
by Gottschalk et al in 1965 for the study of proton-proton bremsstrahlung (Gottschalk 
et al 1965). In this geometry, the two outgoing identical particles are forced to have 
equal energies and polar angles. In our case, when applied to the electron and the 
positron in the final state, the use of the Gottschalk symmetry would mean that k is 
forced to be perpendicular to the plane formed by v and K. This geometry is clearly 
equivalent to intersecting the FDCS with a 2D surface, so as to make it dependent on 
only two parameters. 

Another standard geometry that accomphshes the same reduction is the "energy 
sharing" or coUinear arrangement, where the electron and the positron are forced to move 
along the same direction, i.e. fc_ • fc+ = k^k+. Contrary to Gottschalk's Symmetric 
Geometry, here k is parallel (or antiparallel) to K. Note, that the variation of fe_ and 
fc+ is not independent, since they have to share the energy according to k'^+k'^ ~ v^ — 2e. 
In this paper we will employ this latter geometry. 

4. Vortices in positron-atom ionisation collisions 

Vortices are ubiquitous in classical and quantum systems with large numbers of particles. 
Here, however, we are interested in a different kind of vortices that occur in the wave 
function of a few-body quantum system when considered in the framework of de Broglie 
- Bohm theory (Bohm 1952, Diirr and Teufel 2009) or Madelungs's Hydrodynamic 
interpretation (Madelung 1926, Ghosh and Deb 1982). In 1926 Madelung demonstrated 
that the time evolution of a wave function ^, as given by the Schrodinger equation, is 
equivalent to that of a perfect fiuid of density p = l^&p and velocity u = Im(Vln^) 
described by Euler equations with a pressure tensor 

Pij = V(8 Vlnp 
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The flow is shown to be irrotational except at vortex singularities where the velocity 
diverges. 

In the single ionisation of an atom by positron impact, vortices might be formed in 
the three-body time-dependent wave function \E'(r_|_, r_, t) at the coalescence region i.e. 
when all the intervening particles are close together. They should emerge either as a 
closed submanifold or in pairs of opposite circulation (Bialynicki-Birula et al 2000) and 
collapse at later times, but some of them can eventually survive up to the asymptotic 
regime at macroscopic distances (Macek 2009). In this case, according to the imaging 
theorem (Bollard 1971, Macek et al 2010) 

\Tif \ (X lim t^ \'^(k+t,k_t,t)\ 

they will directly manifest themselves as zeros of the ionisation matrix element Tj/. 

Since both the real and imaginary parts of Tjj have to vanish at a vortex, such 
structure will eventually span a submanifold with codimension 2. By employing a 
collinear geometry, we are reducing the number of independent scalar variables from 
four to two, and therefore vortices should have to show up as isolated zeros. With 
this in mind, we conduct a systematic search of vortices in positron-impact ionisation 
collisions at intermediate energies, by looking for crossings of the projection onto the 
collinear geometry of the sub manifolds where the real and the imaginary parts of Tif 
vanish. Within this geometry, the real and imaginary parts of the normalization factors 
N{Z/kJ) and N{l/2k) vanish on an infinite number of circles accumulating at A;_ = 
and k- = k+. These conditions correspond to the well-known mechanisms of electron 
excitation to the continuum (EEC) and electron-capture-to-the-continuum (ECC), 
respectively. This infinite number of non-crossing lines might mask the emergence 
of a vortex. Thus, in figure [T] we choose to show only the lines where the real and 
imaginary parts of the reduced transition matrix element Tif = Tif/N{Z/k-)N{l/2k) 
vanishes. Obviously, these lines do not coincide with those corresponding to T,/, but 
their crossings do. Thus, their use is instrumental for our task of locating single zeros 
of the ionisation matrix element. As it is shown in figure [T| in the case of the ionisation 
of atomic Hydrogen by positrons of 100 and 200 eV, this only occurs at a single point, 
when the electron is in the close vicinity of the saddle point potential of the positron - 
ion continuum system. 

In figure [T| at the point where the zeros of the real and imaginary part of the reduced 
transition-matrix cross each other, the amplitude vanish. This is clearly observed in 
figure [2| confirming our expectations that the deep minimum previously found by us at 
100 eV (Delia Picca et al 2005) actually represents a zero of the FDCS. Moreover, the 
present results confirm that this is not an isolated effect at 100 eV but that the same 
features are obtained at larger energies. Although some differences in the FDCS are 
observed for the two energies investigated they, present a similar qualitative behavior, 
showing vortices at about the same electron emission direction. 

The generalized velocity field 

u = Im — ^ 
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Figure 1. (color online) Lines of zero real and imaginary parts of the reduced 
transition matrix element Tif (see text) in a coUinear geometry for the ionisation 
of an Hydrogen atom by the impact of a positron of 100 (Left) and 200 eV (Right). fc|| 
and fc_L are the components of the electron momentum fee parallel and perpendicular 
to the initial velocity of the positron, respectively. The axis are normalized to the 
incident momentum fcmax- 
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Figure 2. (color online) Square modulus of the amplitude of the reduced transition 
matrix element (see text) for the ionisation of an Hydrogen atom by the 

impact of a positron of 100 (Left) and 200 eV (Right). Conditions are set to the 
collinear geometry configuration discussed above. fc|| and k± are the components of 
the electron momentum fcg parallel and perpendicular to the initial velocity of the 
positron, respectively. 
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Figure 3. Generalized velocity field of the reduced transition matrix element fif (see 
text) in a collinear geometry for the ionisation of an Hydrogen atom by the impact of 
a positron of 100 eV. 



can be shown to be normal to the submanifolds where the real or the imaginary parts 
of the wave-function vanishes. Thus it produces a circulation along any closed contour 
encircling the vortex that must be quantized to a multiple of 27r in order to assure the 
single valuation of the wave-function (Bialynicki-Birula et al 2000). The velocity field 
associated to the reduced ionisation matrix element Tij is shown in figure [s] In fact, 
the numerical integration of the generalized velocity field around the vortex yielded a 
circulation equal to 2n. 

Further calculations to be shown in a future more comprehensive article let us 
confirm the appearance of other vortices at higher energies, as the one reported by 
Brauner and Briggs (1991). But only the one discussed here remains at intermediate 
and low energies, being visible at energies as low as 30 eV. Furthermore, its location in 
momentum space does not significantly depend on the impact energy. 

5. Conclusions 

In this article we have uncovered the presence of a vortex in positron-impact ionisation 
collisions at intermediate energies that could be observed with reaction-microscope 
or coincident electron-positron spectroscopy techniques. Vortices are ubiquitous in 
classical and quantum systems with large (macroscopic) number of particles. They are 
routinely observed in connection with superfluidity, superconductivity, or Bose-Einstein 
condensation, and are described by introducing nonlinear terms in the Hamiltonian. 
The main qualitative differences of these examples with the vortex investigated in 
this work are, on one hand the extreme simplicity of the present three-body system, 
consisting of an electron, a positron, and a proton. Additionally, the time-evolution 
is governed exactly by the Schrodinger equation with Coulomb interactions, with no 
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ad-hoc potential models or nonlinear terms. In fact, our theoretical prediction might 
prompt the experimental search of positron-induced vortices, not only providing a 
possible confirmation of our findings, but also motivating the study of the topology 
and dynamics of these intriguing quantum structures. 

The vortex investigated here is located in a configuration where the electron sits 
in the vicinity of the positron - ion potential saddle point. In this sense, it is akin 
to a deep minima theoretically observed by Brauner and Briggs (1991) at 45 degrees 
in a similar coUinear geometry with significantly larger impact energies of 10 and 100 
keV. In both cases, vortices seem to be located at significant points related to particular 
collisions mechanisms, as the saddle-point and Thomas processes, respectively. Whether 
this represents a rule or a fortuitous coincidence remains to be investigated. 

Finally, let us point out that all the available studies of vortices in ionisation 
processes, and the one undertaken here is not an exception, have relied on some kind 
of geometry that projects the FDCS onto a two-dimensional plane. This means that 
the isolated zeros observed so far in electron and positron-impact ionisation coUisions 
are only the cut of minimum dimensionality of more complex structure spanning a 
submanifold of codimension 2 in the 4D space of independent scalar variables of 
the FDCS. The future research of this complex structure would clearly represent a 
formidable but inescapable task in order to achieve a fundamental understanding of 
these few-body fascinating structures. 
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